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ABSTRACT 


Jhis  report  consists  of  two  major  parts.  Each  part  addresses 
different  problems  or  issues.  Part  I  is  titled,  "Comparison  of 
"on-board"  In-situ  Vertical  Sound  Speed  Profiles  from  the  SIMAS  and 
I  CAPS  Environmental  Data  Bases  with  High  Quality  CTD  and  STD  Data", 

APP  FY-80  Task  #1.  The  purpose  of  Part  I  was  to  compare  the  SIMAS 
and  ICAPS  sound  speed  profiles  with  sound  speed  profiles  from  assumed 
oceanography  (CTD  or  STD  observations). 

Part  H  is  titled,  "Evaluation  of  the  SIMAS  and  ICAPS  Environmental 
Data  Bases,  Data  Handling  Procedures,  and  Merge  Methodology",  APP  FY-80 
Task  #11.  The  purpose  of  Part  H  was  to  compare  and  recommend  the 
data  bases,  algorithms,  merge  techniques  of  either  SIMAS  or  ICAPS.  The 
major  causes  for  the  "significant"  differences  in  the  SIMAS  vs.  ICAPS 
sound  speed  comparisons  were  found  to  be  in  the  merge  methodologies  and 
historical  environmental  profiles.^. 
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RECOMMENDATIONS 


Based  on  the  findings  of  Part  I  and  II  of  this  report,  the  following 
recommendations  are  forwarded: 

•  Philosophy 

-  Recommend  that  the.  I  CARS  philosophy  of  reliance  on  the  in-situ 
ST  profile  be  adopted.  In  SIMAS,  there  appears  to  be  far  too 
much  reliance  on  the  integrity  of,  its  fixed,  historical  sound 
speed  profiles.  Thu  heavy  reliance  on  a  fixed  so  and  a  peed 
profile  and  a  fixed  cAiteAia  can  lead  to  problems  when  operating 
is  known  physically  dunamic  ocean  regions  where  the  z ones  them¬ 
selves  meander  and  migAate  horizontal,  as  well  as,  vertical 
distances . 

•  In-situ  BT  data  handlina 

-  Recommend  that  the  7CAPS  pAoceduAe  [inputting  of  onlu  observed 
pains  of  depth- temperature  pains)  be  adopted  as  an  "intenim" 
pAoceduAe.  The  SIMAS  pAoceduAes  of  inputting  7-9  inflection 
points  (G.  SAOwn)  by  use  of  a  straight- edge  to  line  up  data 
points,  and  then  picking  off  and  using  as  input  the  numerical 
values  whe ne  the  straight  lines  inteAsect  leads  to  inputting 
depth  and  tempenatune  values  which  may  have  neveA  been  A ecorded. 

-  Recommend  that  A VP  favoAably  encouAage  any  ongoing  ef fonts  to 
inputting  not  selected  data  points  fAom  a  in-situ  BT  pAofile, 
but  the  inputting  of  the  in-situ  BT  tAace  thAough  a  continuous^ 
AecoAding  procedure.  This  will  eliminate  the  variable- inconsistent 
human  facto A  of  incoAAect  selections  of  inflection  points  along 

a  BT  tAace. 

9  Historical  Environmental  Data  Base 

-  Recommend  that  the  ICAPS  historical  tempeAatuAe  and  salinity  data 
base  be  adopted  on  an  interim  basis. 

-  Recommend  two  alternatives  to  replace  or  modify  the  present  TCAPS 
historical  data  bases: 


la)  When  available.,  replace  the.  ICAPS  temperature,  salinity 
{iles  with  the  temperature,  salinity,  density  and  sound 
velocity  analyzed  data  sets  f,rom  the  Generalized  Digital 
Environmental  Model  [GVEM]  by  Dr.  T.  Davis,  (NOG- 022} . 

(b)  Using  recent  deep  X8T  data  las  available)  update  the  archival 
temperature  6 He ,  particularly  in  the  region  o{  the  deep 
sound  channel. 

•  Develop  a  more  appropriate  procedure  to  insure  that  the  proper 
depth  of  the  layer  depths  are  properly  specified  (i.e.,  replacing 
the  present  procedure  of  manual  interpretation  and  inputting  of 

the  in-sltu  temperature  data  with  an  automated  procedure  which  utilizes 
and  inputs  the  entire  trace  or  profile  taken). 

•  Recommend  a  more  appropriate  representation  of  the  depth  of  the  deep 
sound  channel  axis  be  developed. 

•  Additional  Recommendations 

-  Recommend  that  the  bathymetry  {lies  in  TCAPS  be  replaced  by  the 
automated  SW/8APS  ( Synthetic  Bathymetric  Profiling  System) . 
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PART  I 


"Comparisons  of  Vertical  Sound  Speed  Profiles 
predicted  from  the  SIMAS  and  ICAPS  Environmental 
Data  Bases  with  High  Quality  CTD  and  STD  data". 


I.  INTRODUCTION 


The  SIMAS  (Sonar  In-situ  Mode  Assessment  System)  and  ICAPS  (Integrated  Command 
Anti-submarine  warfare  Prediction  System)  programs  address  different  system  appli¬ 
cations  but  are  similar  in  that  each  has  been  developed  for  use  as  an  "on-board" 
prediction  system.  It  is  known  that  SIMAS  and  ICAPS  do  not  employ  identical 
environmental  data  bases.  The  extent  to  which  these  data  bases  differ  is  being 
addressed  in  a  separate  study.  The  purpose  of  this  particular  study  was  to: 

•  compare  predicted  SIMAS  and  ICAPS  vertical  sound  speed  profiles 
with  observed,  measured  and  recorded  in-6-Uu  CTD  (Continuous- 
Temperature-Depth)  or  STD  (Salinity-Temperature-Depth)  data. 

•  present,  whenever  possible,  similarities  or  dissimilarities  between 
the  predicted  SIMAS  and  the  ICAPS  sound  speed  profiles  with  the 
CTD  and  STD  data. 

•  discuss  the  implications  of  these  similarities  and  dissimilarities. 

•  compare  the  "goodness"  of  SIMAS  and  ICAPS  to  the  "Assumed  Oceanography". 

Since  this  study  was  of  limited  scope  it  focused  entirely  on  vertical  sound 

speed  profile  comparisons.  Reliable,  processed  CTD  or  STD  data  available  to 
study  was  (17)  records.  Figure  1  presents  the  geographical  locations  of  the 
test  cases. 

In  an  earlier  study  entitled,  "Comparison  of  the  ICAPS  and  SIMAS  Historical 
Environmental  Data  Bases"  by  E.  Hashimoto,  NORDA  321,  FY-79  APP  Task  II  Report, 
dated  28  September  1979,  the  author  compared  vertical  sound  speed  profiles  from 

v 

the  SIMAS  and  ICAPS  environmental  data  bases  using  XBT  (expendable  bathythermograph) 
data  as  input.  It  also  presented  "significant"  differences  found  between  the  two 
environmental  data  bases.  That  report,  which  described  the  selection  procedures 
of  appropriate  SIMAS  and  ICAPS  historical  sound  speed  profile,  is  being  published 
as  a  NORDA  Technical  Memo  #66  titled,  "Comparisons  of  Historical  and  in-iitu 
Vertical  Sound  Speed  Profiles  from  the  SIMAS  and  ICAPS  Environmental  Data  Bases". 
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In  the  present  study,  SIMAS  and  ICAPS  simulations  of  vertical  sound  speed 
profiles  using  the  STD  temperature  data  are  compared  to  the  measured  sound  speed 
profiles  or  profiles  calculated  by  Wilson's  equation. 


Footnote:  Comparison  of  bottom  loss  provinces  and  bottom  loss  versus  grazing 
angle  curves  used  by  SIMAS  and  ICAPS  were  also  conducted  by  the  author.  At  18 
sites,  ICAPS  and  SIMAS  results  were  compared;  bottom  loss  provinces  were  compared 
In  the  Mediterranean  Sea.  The  results  of  this  study  are  contained  in  NORDA 
Technical  Note  No.  71. 


2 


II.  APPROACH  USED  IN  STUDY 


The  approach  used  In  this  comparison  study  was  to  obtain  and  establish  a 
set  of  test  case  profiles  using  measured  ocean  environmental  data  with  which  the 
predicted  SIMAS  and  ICAPS  sound  speed  profiles  could  both  be  compared.  The  test 
case  profile  data  are  from  reliable  CTD  and  STD  observations.  The  CTD  and  STD 
observations  are  considered  in  this  study  as  being  the  "assumed"  or  real  oceanography. 
CTD  and  STD  test  cases  were  selected  so  as  to  avoid  biasing  the  results  towards 
any  given  system.  The  selections  were  made  on  the  reliability  (author's  confidence 
in  the  observations)  and  availability  (accessible  in  appropriate  format  to  impact 
this  study)  of  CTD  and  STD  data. 

Both  CTD  and  STD  probes  have  an  accurate  temperature  sensor.  The  STD  also 
measures  salinity.  The  values  of  temperature  from  both  sensors  are  recorded 
during  the  observation  period.  The  literature  generally  states  that  the  relative 
accuracy  of  the  CTD  Is  approximately  +  0.005°C  for  temperature.  The  STD  has  a 
relative  temperature  accuracy  of  +  0.02°C,  a  relative  salinity  accuracy  of  +0.02 
parts  per  thousand. 

Identical  measured  values  of  temperature  from  either  the  CTD  or  STD  were  used 
as  input  data  for  both  SIMAS  and  ICAPS  sound  speed  calculation  programs  and 
treated  as  <f  the  values  were  those  taken  by  <n-6ltu  XBTs.  Vertical  sound  speed 
profiles  were  obtained  from  both  systems.  The  SIMAS  and  ICAPS  generated  profiles 
were  plotted  together  with  the  values  of  sound  speed  either  measured  (by  the  CTD), 
or  calculated  from  Wilson's  equation  (for  the  STD).  The  intent  for  such  plots 
was  to  illustrate  the  similarity  between  the  SIMAS  or  ICAPS  vertical  sound  speed 
profiles  and  the  "assumed  oceanography". 

The  designations  used  and  the  locations  of  the  in-i-Uu  test  cases  are  tabulated 
in  Table  I.  The  sources  for  the  CTD  and  STD  test  cases  are  identified  in  Table  II. 
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III.  COMPARISONS 

In  this  section  the  results  of  the  various  comparisons  are  presented  along  with 
supporting  information.  In  the  simulated  comparisons  of  predicte'  SIMAS  and  ICAPS 
produced  sound  speed  profiles  with  "assumed"  oceanography,  similarities  and 
dissimilarities  were  found  with  the  CTD  and  STD  results  (termed  "assumed  oceanography" 
in  table  XV). 

In  some  cases,  similarities  and  dissimilarities  appear  to  be,  more  or  less, 
consistent  or  reoccurring.  In  other  cases,  they  appear  unique  in  their  occurence. 
Presented,  are  the  generalized  similarities  and  major  dissimilarities  found  when 
comparing  the  predicted  SIMAS  and  ICAPS  profiles  with  the  assumed  oceanography. 

The  similarities  and  dissimilarities  are  indicated  in  the  appropriate  sections 
to  follow. 

•  All  temperature  and  salinity  data  from  CTD  and  STD  units  were  sent 

to  NUSC  and  NAVOCEANO,  for  SIMAS  and  ICAPS  processing,  respectively. 

Naval  Underwater  Systems  Center  (NUSC) 

Naval  Oceanographic  Office  (N00) 

•  Results  were  returned  to  NORDA  where  they  were  plotted. 

•  The  accuracy  and  correctness  of  all  plots  were  triple  checked  to 
eliminate  possible  errors  which  may  have  resulted  from  (interpreting) 
reading  the  data  sheets,  key  punching,  unit  conversion,  or  from  plotting 
of  incorrect  data  files. 

•  All  profiles  were  generated  by  personnel  at  NUSC-New  London,  CT  (for  SIMAS) 
and  NOO-Bay  St.  Louis,  MS  (for  ICAPS). 

•  All  "on-board"  input  data  into  SIMAS  and  ICAPS  were  identical 
temperature  profiles  recorded  by  the  CTD  and  STD  temperature  sensors. 

•  Comparisons  for  similarity  and/or  dissimilarities  in  sound  speeds  (! SS ) 

were  conducted  for  values: 

SS,  at  the  surface  (Table  III) 

SS,  at  the  layer  (SSL)  (Table  IV) 

SS,  at  1000  feet  (SS1000)  (Table  V) 

SS,  at  the  sound  channel  axis  (SSCA)  (Table  VI) 
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•  Comparisons  for  similarity  and  or  dissimilarities  in  depth  (D)  were 
conducted  for  values: 

D,  at  the  layer  (in  feet)  (DLD)  (Table  VII) 

D,  at  the  depth  of  the  sound  channel  axis  (in  feet)  (DSCA)  (Table  VIII) 
Figures  2-18  contain  comparisons  of  vertical  sound  speed  profiles  which 
were  generated  from  the  SIMAS  (dashed  lines),  and  ICAPS  (dotted  lines)  environmental 
data  bases,  and  the  in-situ  CTD/STD  (both  as  solid  lines)  observations. 

Tables  (IX),  (X),  (XI),  (XII),  (XIII),  and  (XV)  present  qualitative  differences 
in  sound  speeds  at  the  surface,  at  the  layer,  at  1000  feet,  at  the  channel  axis, 
in  the  depth  of  the  layer  and  the  channel  axis,  respectively.  Table  (XV)  contains 
the  numerical  values  of  sound  speeds  computed  from  CTD  and  STD  sensors.  Table  (XVI) 
contains  the  numerical  values  of  sound  speed  computed  from  SIMAS.  Table  (XVII) 
contains  the  numerical  values  of  sound  speed  computed  from  ICAPS. 
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.  SUMMARY  OF  THE  COMPARISONS 


As  stated  in  the  introduction,  the  purpose  of  this  study  was  to  compare 
predicted  SIMAS  and  ICAPS  sound  speed  profiles  with  (actual)  CTD  or  STD  data 
(considered  as  the  "assumed  oceanography"),  to  present  thr  _imilarities  or 
dissimilarities,  and  to  discuss  any  implication  of  such  similarities  or  dissimi¬ 
larities  for  certain  test  cases. 

Tabulations  of  the  numerical  qualitative  term  and  quantitative  terms  for  SS 
at  the  surface,  at  the  layer,  at  1000  feet,  at  the  deep  sound  channel  axis,  at  the 
depth  of  the  layer,  and  at  the  deep  sound  channel  axis  were  presented  on  tables 
(III  thru  VIII).  There  were  certain  similarities  and  dissimilarities  which  appear 
more  frequently  than  others.  Tables  (IX  thru  XIV)  are  presented  in  an  attempt  to 
simplify  these  similarities  and  dissimilarities.  The  author,  for  SIMAS  and  ICAPS 
and  the  "assumed  oceanography",  has  made  comparisons  of  the  values  in  sound  speed 
(in  ft/sec)  at  the  surface,  depth  of  the  layer,  at  1000  feet,  and  at  the  deep  sound 
channel  axis.  Comparisons  were  made  between  (a)  ICAPS  and  SIMAS  predicted  results  with  (b) 
measured  sound  speeds  from  CTO  or  sound  speed  calculations  from  Wilson's  equation 
from  STD  measurements.  A  brief  summary  of  these  tables  (IX)  through  (XIV)  are  as 
follows: 

For  Sound  Speed  Profile  Compari: on  In  Sound  Speeds 

•  at  the  surface  (17  test  case  comparisons): 

"none  to  slight"  differences: 

ICAPS  =  14  cases 
SIMAS  =  2  cases 

"I  f/s  <SS  <7  f/s"  differences: 

ICAPS  =  3  cases 
SIMAS  =  7  cases 

"7  f/s  <SS*  <20  f/s"  differences: 

ICAPS  =  0  cases 
SIMAS  =  6  cases 

"20  f/s  <SS**  <40  f/s" 

ICAPS  =  0  cases 
SIMAS  =  2  cases 

•  at  the  layer,  (SSL)  (14  test  case  comparisons) 

"none  to  slight"  differences: 

ICAPS  =  14  cases 
SIMAS  =  2  cases 

"1  f/s  <SSL  <7f/s"  differences: 

ICAPS  =  3  cases 
SIMAS  =  7  cases 
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"7  f/s<  SSL*<  go  f/s"  differpnces: 

ICAPS  =  0  cases 

SIMAS  =  6  cases 

"20  f /s <  SSl**<  40  f/s"  differences : 

ICAPS  =  0  cases 

SIMAS  =  2  cases 

•  in  sound  speeds  at  the  channel  axis  out  of  17  test  cases 

"none  to  slight"  di-P-f^rences : 

ICAPS  =  3  cases 

SIMAS  =  0  cases 

"1  f/s<  SSCA<  3  f/s"  differences : 

IlaPS  =  4  cases 

SIMAS  =  2  cases 

"3  f/s<SSCA*<  7  f/s"  differences: 

ICAPS  =  2  cases  ~ 

SIMAS  =  5  cases 

"7  f/s<  SSCA**<  12  f/s"differences: 

ICAPS  =  5  cases  ~ 

SIMAS  =  4  cases 

"12  f/s<  SSCA***<  20  f/s" differences: 

- ICAPS  =  0 "cases 

SIMAS  =  0  cases 

"20  f/s <  SSCA****  "differences: 

ICAPS  =  0  cases  ~~ 

SIMAS  =  3  cases 

•  in  sound  speeds  at  1000  feet  out  of  17  test  cases: 

"none  to  sliqht"  differences: 

— I  CAPS  =  1  {Teases 
SIMAS  =  0  cases 

"1  f/s <  SSI 000 <  6  f/s"  differences: 
it APS  =  5  cases 

SIMAS  =  7  cases 

"6  f/s <  SS1000*''  12  f/s"  differences: 
it  APS  =  2  cases 

SIMAS  =  3  cases 

"12  f/s r  SSI 000** <  20  f/s"  differences: 

- rtAP5~=nr"cases 

SIMAS  =  4  cases 

”20  f/s <  $S1000***<  75  f/s"  differences: 

- ICAPS'  (Teases 

S  If ’AS  =  0  cases 

"75  f/s<  SSI 000**** <  100  f/s"  differences: 

~  'ICAPS' =  O’cases 
SIMAS  =  2  cases 

"100  f/s <  SS.1000*****"  differences: 

1  CAPS  =  0  cases 
SIMAS  =  1  case 
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•  in  the  depth  of  the  layer  out  of  17  cases: 

"none  to  sljohf"  differences: 
rC-fRT  =  3  cases 
SIMAS  =  3  cases 

"2  feet v  OLD < 30  feet"  differences: 
rCAP5~=  9  cases 
SIMAS  =  8  cases 

“30  feof<  DL3*<  90  feet"  differences: 

ICAPS  =  3  cases 
SIMAS  =  4  cases 

"90  feet<  OLD**  <150  feet"  differences: 

I  CAPS  =  1  case 
SIMAS  =  0  cases 

"150  feet<  DLD***<250  feet"  differences: 
ICaPs  =’l  case  ~ 

SIMAS  =  0  cases 

"250  feet<  DLC****"  differences: 

~ ~  I CAPS  =“0'ca s e s 
SIMAS  =  2  cases 


•  in  the  depth  of  the  primary  sound  channel  axis,  out  of  14  test  cases 

"less  than  300  feet"  differences: 

ICAFS~=~7“ cases 
SIMAS  =  6  cases 

"3QC  feet  <  DSCA* <  GOG  feet"  differences: 

ICAPS  =  3  cases 
SIMAS  =  2  cases 

"600  feet < DSCA** <  1200  feet"  differences: 

ICAPS  =  4  cases 
SIMAS  =  2  cases 

"1200  feet  < DSCA***  < 2000  feet"  differences: 

ICAPS  =  0  cases 
SIMAS  =  3  cases 

"2000  feet  < OSCA****  "  differences: 

ICAPS  =  0  cases 
SIMAS  =  1  case 


•  ICAPS  utilized  the  CTD/STD  input  temperature  data  above  1000  feet 
more  than  SIMAS. 

•  The  "rejection"  of  in-iiXa  temperature  information  in  SIMAS  is  clearly 
too  restrictive.  The  consequent  use  of  archival  information  results  in 
large  discrepancies  with  CTD/STD  results. 
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VI .  CONCLUSIONS 


In  t**e  previous  study  (NORDA  TM  No.  66)  which  compared  the  SIMAS  and  ICAPS 
Historical  and  IWSITU  Environmental  Data  Bases  vertical  sound  speed  profiles  and 
presented  "significant"  differences,  there  was  no  comparisons  made  with  real  IWS7TU 
measurements.  In  this  study,  this  has  been  conducted  to  test  for  "goodness"  and 
the  results  are  presented. 

Upon  review  of  the  comparison  plots  and  the  tables  of  similarities  and 
dissimilarities,  it  is  clear  that  ICAPS  demonstrates  a  noteable  superiority  to 
SIMAS  in  its  ability  to  predict  measured  sound  speed  from  an  "XBT"  obtained  from 
CTD  or  STD  data.  This  does  not  imply,  however,  that  ICAPS  does  not  need  improvement 
as  will  be  discussed  under  recommendations. 

The  authors  conclusions  based  on  tables  (III  through  VIII)  and  (IX  through  XIV) 
are  as  follows: 

•  "Assumed  oceanography"  was  reproduced  more  closer  by  ICAPS. 

t  Sound  speeds  at  the  surface  were  reproduced  more  closer  in  the  North 
Atlantic  and  North  T a  ' f i c  Oceans  by  ICAPS. 

ICAPS  AefZectt  a  tZgnifieantZy  ZaAgeA  quantity  (14  out  of,  17}  of, 
caaea  wheAe  tl le  vaZuet  o f  tound  apeeda  at  the  tuAface  ag4.ee  [<  1  i(/a) 
cZoteZy  with  the  CTV/STV  obteAvationt  than  SIMAS  (2  out  of  17). 

Whene  ICAPS  hat,  zeAo  caaea,  SIMAS  thowed  8  out  of  17  coaea  wke Ae 
the  vaZuet  of  tound  apeed  at  the  tuAface  diffeA  between  eitheA  7  to  20  i(/a, 
oa  20  to  40  i (/a. 

•  Sound  speeds  at  the  layer  were  reproduced  more  closer  in  the  North  Atlantic, 
North  Pacific  and  Indian  Oceans  by  ICAPS. 

ICAPS  AefZ.ectt  a  tignificantZy  laAgeA  quantity  [14  out  o f  17)  of 
caaea  whe Ae  the  vaZuet  of  tound  apeeda  at  the  ZayeA  ago.ee  (< 1  f/t)  cZoteZy 
with  the  CTV/STV  obteAvationt  than  SIMAS  (2  out  of  17). 

WheAe  ICAPS  hat  zeAo  caaea,  SIMAS  thawed  8  out  of  17  catet  wheAe 
the  vaZuet  of  tound  tpeed  at  the  ZayeA  diffeA  between  eitheA  7  to  20  f/t, 
oa  20  to  40  f/t.  10 


•  Sound  speeds  at  1000  feet  were  reproduced  more  closer  in  the  North 
Atlantic  Ocean  by  ICAPS. 

ICAPS  fie.ile.cti  a  iigrUficantly  laAgeA  quantity  ( 10  cut  of  17) 
of  ana  whene  the  value*  of  iound  ipeedi  at  1000  feet 

agnee  (<lf/i)  cloieA  with  the  CTD/STV  obiefivationi  than  S7MAS 
(0  oat  of,  17).  SIMAS  ihowed  7  caiei  to  be  g neaten  than  12  f/i 
at  1000  feet  of  which  3  of  thii  7  afie  gneate A  than  75  f/i.  ICAPS 
on  the  otheA  hand  ihowi  all  of  iti  value*  at  1000  feet  not  exceeding 

n  ft 6. 

•  Neither  ICAPS  nor  SIMAS  were  impressive  in  their  handling  of  sound  speeds 
at  the  sound  channel  axis. 

ICAPS  ihowed  that  In  9  out  of  17  caiei,  the  value i  of  iound 
ipecd  at  the  axii  wene  *7  f/i  and  SIMAS  ihowed  7  out  of  17  wene 
•<7  f/i.  SIMAS  ihowed  3  caiei  to  be  >20  f/i  at  the  channel  axti, 
and  ICAPS  had  none. 

•  Neither  ICAPS  nor  SIMAS  were  impressive  in  their  handling  of  the  depth  of 
the  layer. 

ICAPS  ihowed  the  depth  of  the  latjem  to  diffen  between  7  to  30  feet 
in  9  out  of  17  caiei  and  SIMAS  ihowi  the  depth  of  the  lay cm  to  differ 
betu'cen  7  to  30  feet  In  S  out  of  17  caici. 

•  Neither  ICAPS  nor  SIMAS  were  impressive  in  their  handling  of  the  depth  of 
the  sound  channel  axis. 

ICAPS  ihowed  the  depth  of  the  iound  channel  axii  to  diffen. 
between  0  -  3 00  feet  in  7  out  of  14  caiei,  and  SIMAS  ihowed  it  in  6  out 
of  14  caiei.  ICAPS  ihowed  the  depth  of  the  iound  channel  axii  to 
diffen.  by  gneateA  than  600  feet  in  4  out  of  14  caiei  and  SIMAS  ihowed 
it  in  6  out  of  14  caiei. 
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GEOGRAPHICAL  LOCATIONS  OF  THE  SITES  FOR  COMPARISON 


pth  (in 


0000 


AO -05 


BOTS 


000 8 


B0Ml  00001  0000 


0000 


in  Ft) 


0000 


in  Ft) 


sm 


Sound  Speed  (in  Ft/Sec) 


000S 


2000  5000  0000  8000  10000  1 2000  1 5000  1  6000 


Sound  Speed  (in  Ft/Sec) 


eiOd 


4000  8000  8000  10000  12000  M000  10000 


TABLE  (I)  DESIGNATORS  &  LOCATIONS  OF  SELECTED  SITES 

FOR  SIMAS,  I  CAPS  AND  CTD/STD  COMPARISONS 


North  Padific  Ocean 


. 

Latitude 

Longi tude 

Month 

P01 

43d  00 ‘ N 

158P02*W 

(September) 

P02 

35°59'N 

•r 

158°0Q'W 

(October) 

P03 

27°20'N 

154°01  ’  E 

(Apri 1 ) 

P04 

33°18’N 

I54°03’E 

(Apri 1 ) 

P05 

36°40'N 

1 54°00 ' E 

(April ) 

P06 

35°00‘N 

I 30°1 .  4 '  W 

(March) 

P07 

35°00’N 

137°59'W 

(March) 

P08 

34°56 . 5 1 N 

1 47°59. 4 1 W 

(March) 

P09 

35°I 2 . 2 1 N 

147°0.3'E 

(April ) 

PI  0 

34°58.2'N 

1 46°01 ' E 

(April ) 

Indian  Ocean 

101 

4°49'N 

053°02'E 

(March) 

N orth  Atlantic  Ocean 


A01 

37°00'N 

01 8°00'W 

(February) 

A02 

36°00'N 

019°00’W 

(February) 

A03 

35°22 . 7 1 N 

035°4 . 7 '  W 

(June) 

A04 

35°12.4'N 

054°59. O' W 

(July) 

A05 

30°24 . 6 ' N 

071°57 . 8' W 

(September) 

A06 

30°25.6'N 

071°57. O' W 

(September) 
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TABLE  (II) 

Sources  for  the  CTD  and  STD  Test  Cases 


I .  Scr  ipts  Ins t i  jpte  for  _Ocean_op_r_aph^_ 

a)  CTD 

j.  Reid  and  K.  Kenyon  P01 ,  P02,  P03,  P04,  P05 

b)  STD 

J.  Reid  and  K.  Kenyon  P06 ,  P07,  P08,  P09,  PI 0 


II.  Woods  H o 1  e_  0 c ea nographic  Institute 
a)  CTD 

M.  McCartney  A05,  A06 


III.  Naval  Ocean  Research  and_Dpv_eJo_pinen^J\ctivj_t^ 
a)  CTD 

H.  Perkins  and  D.  Fenner  A01,  A02,  101 

IV.  Naval  Oceanoqraphi c  Office 
a)  CTD 

LCDR  L.  Danzler  A03,  A04 


31 


e,  a 


)  at  the  surfac 


Differences  in  ?s  (in  ft /sec) 
at  the  surface 


SIMAS  At '.lined  ’.CAPS -Auuned 


Oceanography 

10.00 


12.03 


-1.47 


-2.31 


-0.80 


-0./8 


-l.il 


12.04 


-4 .  20 

30.45 


-13.42 


-2.66 


5.78 


-10.46 


-6.66 


Oceenouraphy | 

not  applicable ] 


.65 


-0.33 


-0.98 


-0.65 


-0.33 


0.99 


2.30 


0.98 

0.66 


.0.98 


-0.99 


-3.25 


0.14 


•0.  33 


72.84 


0.45 


i  feet/sec)  at  1000  feet,  and  quant 
feet. 


f  t/‘.ec) 

D)  f  in 

at  1000  fppt 

‘S  (  ir.  ft. 

A  buried 
Ocoancqraphy 

s,v, <  jr ed 
Ocr-anoqrjphy 

KM'C  -t'.s 

4936.62 

9.32 

-0.12 

44<5.75 

19.19 

0.74 

4943.60 

16.20 

0.32 

t:y)A  *] 

-2.61 

-0.65 

9000. 1 1 

-1.11 

0.29 

9001.10 

>2.10 

. . . 

0.36 

38 

-5.45 

0.22 

4690.10 

18.85 

2.18 

4969.24 

-3.30 

0.57 

48%.  54 

94.66 

3.10 

4  <i06. 91 

19.96 

?.  94 

4671.22 

-6.29 

-0 . 69 

48P4.96 

4.97 

-0.76 

4963.18 

16.44 

1 .26 

4HR0.69 

89.02 

6.99 

4353.29 

116.44 

2.97 

4929.65 

-1.65 

7.57 
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TABLE  VIII.  Numerical  values  of  the  channels  axis  depth  (in  feet),  quantitative 
differences  in  the  depths  of  the  channel  axis  (in  feet.) 


table  (XI  ) 


LEGEND:  si iqht  <  1  f/s 

1  f/c <  SSIOOO <  6  f/s 
6  f/s  <  SSIOOO’  <  12  f/s 
12  f/s <  SSIOOO*'  <  20  f/s 
20  f/s  <  SSIOOO*'*  <  75  f/s 
75  f/s<  SSIOOO***’  <  TOO  f/s 
100  f/s<  SSIOOO***** 


40 


TABLE  (xii) 


Oualitative  Differences  in  SS 
at  the  channel  axis 


Differences  in  SS  at 
the  channel  axis 


SITE  (ID) 

SIMAS 

ICAPS 

A01 

SSCA** 

SSCA* 

A02 

SSCA** 

SSCA** 

A03 

-- 

1 

AO  4 

-- 

A  05 

SSCA* 

I  SSCA* 

A  06 

,  SSCA* 

SSCA* 

P01 

.  .  - 

P02 

SSCA* 

SSCA** 

P03 

SSCA* 

SSCA 

PQ4 

SSCA**** 

si  ight 

P05 

SSCA* 

SSCA 

P06 

SSCA 

slight 

P07 

SSCA 

SSCA** 

P  u8 

SSCA** 

SSCA** 

P09 

SSCA**** 

SSCA 

P10 

SSCA**** 

slight 

101 

SSCA** 

SSCA** 

LEGEND : 


-  -  r  no  comparison  made 
s  1  i  n  h  t.  1  f  /  s 
1  f/s  cSfA.-  3  f/s 
3  f/s  ''  SSCA-  -  1  f/s 
7  f/<-  ■  ssf A*’  '12  f/s 
12  f/s  ••  S SPA ?0  f/s 
20  f/v  •'  cSCA**** 


TABLE  (XIII) 


Qualitative  Differences  in  SS 
at  the  layer  depth 


SITE  (ID) 


LEGEND:  si iciht <  2  ft 

2  ft <  QLD<  30  ft 
30  ft<  DLD*<  90  ft 
90  ft<  QLD**<  150  ft 
150  ft<  DLD***<  2o0  ft. 
250  ft<  OLD***'* 


Differences  in  layer  depth 
(in  feet) 


SIMAS 

I  CAPS 

[}[_">**** 

DLD 

OLD**** 

Old*** 

si  iciht 

none 

DLD 

DLD* 

OLD 

OLD* 

•  si ight 

none 

DLD 

DLD* 

DLD 

DLD 

none 

none 

DLD 

DLD 

DLD* 

DLD** 

DLD 

DLD 

DLD* 

DLD 

DLD* 

DLD 

DLD 

DLD  1 

DLD* 

DLD 

DLD 

DLD 

TABLE  (XIV) 


Oualitative  Differences  in  the  depth 
of  the  cuanne1  axis 


LEGEND: 


Differences  in  the 

decth  of  sound  channel  axis  (in  feet) 


- =  no  corperisons  n^ade 

0  ft<  DSC A<  “500  ft 

300  ft<  nsc.p*<  600  ft 

600  ft<  DSCA**<  1200  ft 
1200  ft<  DSCA***<  2000  ft 
2000  ft  <  DSCA**** 


SIMAS 


I  CAPS 


DSC  A*  DSC A** 


DSCA* 


DSC  A* 


DSCA** 


DSCA*' 


DSCA***  OSCA 
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TABLE  (XVI).  Numerical  values  of  sound  speeds  (in  feet/sec)  for  SIMAS  profiles 
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I 


PART  (II) 

"EVALUATION  OF  THE  SIMAS  AND  ICAPS 

ENVIRONMENTAL  DATA  BASE,  DATA  HANDLING  PROCEDURES  AND  MERGE  METHODOLOGY" 


* 


I.  INTRODUCTION 

The  SIMAS  (Sonar  ^"Situ  Mode  Assessment  System)  and  ICAPS  (Jntegrated  Command 
ASW  prediction  System)  environmental  data  bases  are  both  being  developed  and 
available  to  APP.  Each  system  has  been  designed  and  developed  by  different 
Navy  activities  and  address  different  system  applications.  The  common  denominator 
between  SIMAS  and  ICAPS  is  that  both  utilize  an  environmental  data  base  of 
oceanographic  parameters,  both  rely  on  an  in-situ  XBT  Expendable  bathythermograph) 
as  input,  and  both  have  been  developed  as  "on-board"  prediction  systems. 

As  requested  by  the  APP  Data  Base  Manager,  NORDA  321  has  been  tasked  to 
identify,  as  much  as  possible  the  reason(s)  for  the  "significant"  differences  in 
the  final  computed  sound  speed  profiles  from  previous  SIMAS  and  ICAPS  comparisons. 
In  this  portion  of  the  report,  the  test  comparisons  and  evaluations  designed  to 
identify  the  specific  cause(s)  or  reason(s)  for  differences  between  the  SIMAS 
ICAPS  computed  sound  speed  profiles  are  discussed. 


1 1 .  BACKGROUND 


In  FY-79,  a  study  ("comparison  of  the  ICAPS  and  SIMAS  Historical  Environmental 
Data  Base",  by  E.  Hashimoto,  dated  28  September  1979),  was  conducted  to  compare  the 
SIMAS  and  ICAPS  computed  sound  speed  profiles  at  selected  geographical  locations, 
for  a  particular  season  or  month.  From  this  point,  that  report  will  herein  be 
referred  to  as  the  NORDA  TM  No.  66.  The  purpose  of  that  study  was  to 

present  and  compare  any  "significant"  differences  found  between  the  SIMAS  and  ICAPS 
final  computed  sound  speed  profiles.  The  evaluation  and  cause  for  differences 
found  specifically  was  not  to  be  addressed  in  that  study. 

Since  "significant"  differences  were  found  in  the  final  computed  vertical 
sound  speed  profiles,  NAVSEA  06H4  and  NORDA  534  (during  an  APP  sponsored  meetinq  on 
21  Feb  1980)  tasked  NORDA  321  to  conduct  follow-on  studies  to  evaluate  the  ability 
for  SIMAS  and  ICAPS  to  reproduce  real  or  assumed  oceanography,  and  identify  the 
cause(s)  for  previous  differences.  Part  I  of  this  report  presented  the  ability 
of  SIMAS  and  ICAPS  to  reproduce  the  assumed  oceanography  by  using  input  temperature 
data  from  very  accurate  CTD  (continuous-temperature-depth)  or  STD  (salinity-temperature- 
depth)  recorders. 

The  identification  of  possible  cause(s)  or  reason(s)  for  previous  differences  found 
between  the  SIMAS  and  ICAPS  computed  vertical  sound  speed  profiles  will  now  be  presented 
and  discussed 
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F  —  1  1  '■  - 

I 

HI.  APPROACH  USED  IN  THIS  STUDY 

In  order  for  the  author  to  properly  determine  or  reasonably  explain  the 
"probable"  cause(s)  associated  with  the  environmentally  "sionificant"  differences 
found  in  the  APP  FY-79  Task  #2  Report,  a  very  large  effort  was  placed  towards 
understanding  several  key  areas: 

t  The  review  of  all  related  documentation  considered  relevant  to  the 
profile  generation  process  for  each  system  (i.e.,  technical  reports, 
software  coding,  etc.), 

t  Personal  conversations  with  Cognizant  individuals  at  the  Naval 
Underwater  Systems  Command  (NUSC  -  E.  Podeszwa,  Code  3351)  and  the 
Naval  Oceanographic  Office  (N00  -  P.  Moresdorf,  Code  9 200). 

0  The  identification  of  the  sources  for  the  oceanographic  data  used  in 
the  development  of  each  historical  base. 

0  The  identication  of  the  analysis  procedure(s)  used  in  the  selection 
of  historical  sound  speed  profiles. 

0  The  identification  of  the  spatial  and  temporal  resolutions  used  in 
each  system. 

0  The  identification  and  influences  of  the  external  and  internal  methodoloay 
employed  by  each  system  regarding  the  handling  of  in-situ  XRT  data 

0  The  identification  and  influence  of  the  selection  criteria  and  retrieval 
procedures  for  the  proper  historical  sound  speed  profile. 

0  The  identification  and  influence  of  the  merging  techniaues  employed  when 
combining  the  in-situ  XBT  with  the  historical  deep  Drofile. 

0  The  identification  and  influence  of  the  different  sound  speed  equations. 

The  author  had  requested  assistance  from  the  Ocean  Data  Systems,  Inc.  Rockville, 
MD  (Mr.  J.  Lock! in  and  Mr.  B.  W.  Scaife).  This  request  was  approved  by  the  APP  Data 
Base  Manager.  ODSI  was  requested  to  provide  supportive  in-depth  investigations  as 
required. 
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IV.  THE  STUDY 

This  section  will  describe  the  specific  work  conducted  and  material  considered 
to  properly  address  the  issues  presented  in  section  III.  In  order  to  assist  the 
reader  in  understanding  the  various  investigations  conducted,  this  section  has  been 
divided  into  several  categories. 

0  differences  in  philosophies 

0  differences  in  Input  BT  data  handling  procedures 
•  differences  in  the  historical  environmental  data  base  composition 
0  differences  in  the  merge  methodoloqy  of  BT  to  historical 
0  differences  in  the  application  of  salinity 
0  differences  in  the  sound  speed  equations 
0  differences  in  the  software  program  coding  errors 

A .  PhUoiopkici 

There  appears  to  be  fundamental  differences  between  I  CAPS  and  SIMAS 
in  their  basic  conceptual  philosophies: 

SIMAS  =  The  entire  SIMAS  philosophy  rests  on  the  integrity  of  the 

historical  sound  speed  profiles  in  the  environmental  data  base. 

ICAPS  =  The  entire  ICAPS  philosophy  rests  on  the  integrity  of  the 
in-iltu  BT. 
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8.  INPUT  BT  VAT  A  HAHVLWG  PROCEDURES 

Both  SIMAS  and  IGAPS  utilize,  as  an  initial  in-situ  input,  a  BT  (bathythermograph) 
trace.  Each  system  attempts  to  capture  the  salient  features  in  the  temperature 
profile  by  inputting  selected  pairs  of  depth  and  temperature  along  the  trace. 

Although  the  philosophy  to  utilize  an  in-situ  BT  observation  are  similar,  each 
system  employs  a  different  subjective  procedure  which  determines  what  pairs 
of  depths  and  temperatures  will  be  used  as  input. 

The  following  are  brief  descriptions  on  the  BT  data  handling  procedures 
employed  by  each  system: 

•  SIMAS 

Certain  "significant"  points  of  inflection  along  each  BT  profile  are 
marked  or  noted  by  the  SIMAS  operator.  These  certain  "significant"  points 
of  inflection  are  selected  according  to  how  well  they  (the  observed  data 
points)  group  themselves  in  a  more  or  less  a  linear  fashion  (personal 
conversation  with  G.  Brown,  NUSC).  The  determination  irr  the  degree  of 
linearity  of  the  data  points  along  a  BT  trace  is  accomplished  by  using  a 
straight  edge  ( i . e. ,  ruler)  to  line  up  any  data  points  which  fall  along 
the  straight  edge.  While  proceeding  down  the  trace,  if  the  linear  distribution 
and  direction  of  the  observed  data  points  should  change,  the  straight  edge 
is  then  reoriented  in  a  new,  different  direction  (personal  conversation  with 
G.  Brown,  NUSC).  An  example  of  this  process  is  illustrated  on  figure  1. 

According  to  the  personnel  who  have  designed  the  SIMAS  system,  they  prefer 
that  the  maximum  number  of  inflection  points  be  held  to  7  points  per  input 
BT  trace.  As  a  result,  all  data  points  describing  a  given  BT  may  not  be 
inputted  into  SIMAS  (e.g.,  If  profile  has  12  inflection  points  and  a  maximum 
of  7  are  selected,  your  significant  points  could  be  omitted). 


•  I CAPS 


Certain  "significant"  points  of  inflection  along  each  BT  profile 
are  marked  or  noted  by  the  ICAPS  operator.  The  first  depth-temperature 
pair  inputted  into  ICAPS  is  the  surface  and  the  surface  temperature. 

The  remaining  pairs  cf  depth-temperature  are  selected  by  the  operator 
according  to  inflection  which  exist  in  each  BT  trace  where  sudden  changes 
in  trend  (positive  to  negative,  negative  to  positive)  or  gradient  (numerical 
value  per  unit  depth)  occur.  The  numerical  values  of  depth-temperature 
must  actually  occur  and  be  contained  on  the  BT  trace.  All  values  inputted 
are  observed  values. 
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C.  Historical  Environmental  Data  Base  Composition 


Both  environmental  data  basts  have  been  developed  using  the.  ocean  station 
data  observations  at  the  National  Oceanographic  data  Center  (NOVC) . 

Each  system  utilizes  an  environmental  data  base  in  two  modes:  with 
an  in-situ  XBT,  or  without  an  in-situ  XBT  [pure  historical) .  Briefly  presented 
below  are  the  parameters  on  kite  and  available  in  STMS  and  7 CAPS: 


ST  HAS 

7  CAPS 

COMMENTS 

Temperature 

NO 

yES 

TCAPS  Temperature  profiles 
are  surface  to  bottom. 

a)  monthly 

NO 

NO 

b )  seasonally  averaged 

NO 

YES 

TCAPS  Temperature  profiles 

c)  annual 

NO 

NO 

are  seasonal  surface  to 
bottom. 

Salinity 

NO 

YES 

S7MAS  uses  a  constant 
salinity  35  ppt  surface 
to  the  bottom. 

a}  monthly 

NO 

NO 

TCAPS  salinity  profiles 

b)  seasonally  averaged 

NO 

YES 

are  seasonal  surface  to 
bottom. 

Density 

NO 

NO 

a}  monthly 

NO 

WO 

b)  seasonally 

NO 

NO 

Sound  speed 

YES 

WO 

TCAPS  computes  all  value i 
of  sound  speed 

a)  monthly 

YES 

NO 

b)  seasonally 

NO 

NO 

STMAS  sound  speed  profiles 

are  surface  to  bottom. 

Bathymetry 

7 n  addition  to  sound  speed,  temperature  and  salinitu,  there  is  another 
parameter  which  is  common  to  both  ST  NAS  and  7  CAPS .  This  parameter  is  batfametru 
(ocean  bottom  depth}.  Brieflu  described  below  are  the  characteristics  and  sources 
far  the  bathymetry  used  in  each  system ; 
o  I CAPS  (HEATON-NOO} 

Automated,  bathymetric  depths  recorded  on  the  data  f lies  are  those 
of  the  deepest  identifiable  bathymetric  contour  interval  found  in  each 
7 / 4- quadrangle  (every  Ml*  x  Ml *  square  or  3 O'  x  30'}. 


Source: 

-  Mo  nth  Atlantic  and  Mediterranean  In  W  aval  Cc.tayiogfw.phia 
O^ice.  Special  Publication  *1304,  1975  [Confidential) 

-  North  Pacific  ii  Naval  Oceanographic  Office  Special 
Publication  * 1301-2-3 ,  1973  ( Confidential ). 

-  North  Indian  Ocean  i*  WPCS  Chant* ,  Second  Edition, 
o  SIMAS  [Pode&zwa  -  NUSC) 

Manual  point  on.  anea  extraction  off  bathymetric  contour  chant*. 
Source: 

-  North  Atlantic  and  Mediterranean  it>  Naval  Oceanographic 
Office  Special  Publication  *1304,  1975  ( Confidential ) 

-  Nofith  Pacific  ii  Naval  Oceanographic  Office  Special 
Publication  * 1301-2-3 ,  1973  [Confidential.) . 

-  North  Indian  Ocean  i*  WPCB  Chart*,  Second  Edition. 
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V.  MERGE  MET HOVOLOGV 


Each  system  computes  sound  speed  as  a  final  parameter.  However,  a  basic 
difference  exists  in  the  manner  in  which  SIMAS  and  ICAPS  each  merge  to  the  In-si tu 
BT  observation.  The  basic  different  is  merging  to  a  sound  speed  profile  or 
merging  to  a  temperature  profile. 

An  outline  description  of  the  differences  in  the  ICAPS  and  SIMAS  methodology 
of  generating  sound  speed  profiles  is  illustrated  below: 

METHODOLOGICAL  FEATURES  OF  SSP  GENERATOR  (OPS I) 

ICAPS  SIMAS 


Multiple  water  mass  (up  to  5)  capability 
for  each  geographic  area. 

Historical  data  base  contains  for  each 
water  mass,  temperature  and  salinity 
profiles  (ZHi,  THi ,  SHi),  i=l ,  NH^45. 

Extend  the  Bathythermograph  (BT)  profile, 
(ZBi,  TB-j ) ,  j  =  l,  NBT,  with 
historical  temperatures  for  corresponding 
water  mass  at  depth  below  the  last  BT 
and  interpolate  salinity  values  from 
historica.  profile. 


A  geographic  area  is  represented  by  one 
water  mass. 

Historical  data  base  contains  for  each 
water  mass  a  sound  speed  profile  (ZH-j, 

VH,),  NH£27,  and  a  representative  near¬ 
surface  salinity  value  ,  S. 

Utilize  Leroy's  equation  to  calculate  the 
SSP  for  BT  profile  with  historical  constant 
salinity.  The  BT  profile  may  consist  of  up 
,  to  24  points. 


Insert  bathymetry  and  truncate  SSP  at 
bottom  depth. 


Utilize  Wilson's  equation  to  calculate 
SSP. 


Shift  the  BT  SSP  by  the  linearly  interpolated 
velocity  difference  between  historical  SSP 
and  BT  SSP  at  1000  ft.  Next,  extend  below 
the  BT  SSP  with  unmodified  historical  SSP 
points. 

Insert  bathymetry  and  truncate  SSP  at 
bottom  depth. 


Calculate  the  layer  depth. 


Adjust  fathometer  bathymetry  and  interpolate 
sound  speed  at  adjusted  bottom  depth. 


Calculate  the  layer  depth. 


Graphical  presentations  of  the  above  outline  description  are  given  in  the  Logic 
Flow  Diagrams  of  Appendix  (  A  )  and  (  B  )  for  SIMAS  and  ICAPS,  respectively. 
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•  ICAFS  =  The  numerical  values  of  temperature  from  the  In-situ  BT  observation 

are  inputted  into  various  sets  of  temperature  criteria  (as  explained 
in  Appendix  (B)).  After  the  in-situ  BT  profile  has  met  certain 
established  temperature  criteria,  it  is  mexgtd  with  an  analyzed 
seasonal  tempenatufie.  profile  contained  in  the  ICAPS  environmental 
data  base  (see  Appendix  G).  This  merge  temperature  profile  is  then, 
inputted  along  with  an  analyzed  seasonal  salinity  profile  into  Wilsc's 
equation  A  derived  sound  speed  profile  is  then  computed.  At  no  time 
is  the  BT  profile  rejected  or  not  used  once  it  is  inputted  into  the 
ICAPS  ‘-ystem. 

•  SIMAS  =  The  numerical  values  of  temperature  from  the  in-situ  BT  observation, 

along  with  a  constant  salinity  profile  (value  of  35.00  parts  per 
thousand),  are  both  entered  into  Leroy's  equation.  A  derived  sound 
speed  profile  is  then  inputted  into  various  sets  of  sound  speed 
criteria  (as  explained  in  Appendix  (A)).  After  the  derived  sound 
speed  profile  has  met  certain  established  sound  speed  criteria,  it 
is  mcAged  with  an  analyzed  monthly  averaged  ioamd  ipecd  profile 
contained  in  the  SIMAS  environmental  data  base.  If  the  derived 
sound  speed  profiles  does  not  meet  the  sound  speed  criteria,  it  is 
possible  for  that  particular  BT  observation  and  the  derived  sound 
speed  profile  to  be  Ae/ec-ted  and  not  used  for  failure  to  meet  the 
criteria. 

Upon  rejection  of  an  insitu  XBT,  a  revised  sound  speed  profile  is 
computed.  This  revised  sound  speed  profile  is  basically  a  modifi¬ 
cation  of  the  SIMAS  historical  sound  speed  profile.  According  to 
investigations  by  Locklin  and  Scaife  (ODSI),  the  historical  SIMAS 
sound  speed  proPile  is  revised  as  follows: 

1.  The  value  of  the  surface  sound  speed  from  the  historical 
profile  is  used. 
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2.  The  value  of  the  sound  speed  at  the  layer  from  the  historical 
profile  Is  used. 

3.  The  value  of  the  depth  of  the  layer  is  revised:  depth  of  the 
layer  becomes  the  value  of  the  second  depth  point  off  of  the 
XBT  trace.  In  other  words,  the  new  depth  of  the  layer  Is  the 
depth  of  the  XBT  layer. 

4.  Values  of  sound  speed  below  the  revised  layer  depth  are  those 
from  the  historical  sound  speed  profile. 

0  cauaed  by  cU^eM.e.nt  meAg-tng  ttchniquu 

As  stated  above  (and  described  in  Appendix  (B))  ICAPS  does  not  reject 
the  In-situ  BT  at  any  time.  The  values  of  temperature  at  the  surface  and 
and  of  temperature  and  depth  at  the  boundary  of  the  surface  layer  are  inputted 
into  ICAPS  and  maintained  unaltered  throughout  ICAPS  processing.  SIMAS 
on  the  other  hand  will  reject  the  entire  BT  if  the  sound  speed  criteria 
is  not  met.  Upon  rejection  of  a  BT,  the  SIMAS  system  generates  a  revised 
sound  speed  profile  for  use.  Below  is  a  comparison  of  the  differences  in 
surface  sound  speeds  and  layer  depths  due  to  the  different  merging  techniques 
when  using  identical  BT  as  input.  The  major  differences,  as  noted  by 
an  asterisk  (*)  are  cases  where  SIMAS  rejected  the  in-situ  BT  data  and 
computed  revised  sound  speed  profiles.  The  sites  (1C,  IF,  IF,  etc.  are 
test  sites  from  Fig.  2  (taken  from  NORDA  TM  #66). 
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r  i 

ICAPS  AND  SIMAS  LAVER  DEPTHS  AND 
SURFACE  SOUND  SPEEDS  USING 

;  IDENTICAL  BT  AND  HISTORICAL  PROFILES 

j  (DIFFERENCES  CAUSED  BY  MERGING) 


Site 

ICAPS 

Layer  Depth 
(ft) 

Surface 

Sound  Speed 
(ft/sec) 

SIMAS 

Layer  Depth 
(ft) 

Surface 
Sound  Speed 
( ft/sec) 

1C 

469.2 

4964 

469.2 

4943* 

IF 

357.6 

4910 

357.6 

4881* 

IF 

65.6 

4936 

65.6 

4944 

2A 

85.3 

5046 

85.3 

5038* 

2E 

3281.0 

4903 

1220.5 

4909 

2E 

72.2 

4987 

72.2 

4982 

2H 

420.0 

4912 

52.5 

4798* 

3A 

59. 1 

5059 

59.1 

5063 

o  The  historical  6 mind  speed  profiles  for  both  SIMAS  and  ICAPS  were  computed  bu 
the  national  Oceanographic  Data  Center  (HCVC]  using  Wilson's  equation.  Although 
ICAPS  elects  to  continue  using  Wilson's  equation  during  the  conversion  the 
in-situ  ST  pro  pies  to  a  sound  speed  profile,  SIMAS  switches  to  Lercif's  equation 
for  its  conversion  to  a  sound  speed  profile. 
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E.  APPLICATION  OF  SAL7N1TV 


Each  system  employs  different  approaches  in  their  application  of  salinity 
to  the  computation  of  sound  speed.  Obvious  differences  exist.  For  example,  SIMAS 
uses  a  constant  value  of  35.00  parts  per  thousand  (surface  to  the  bottom)  and  I  CAPS 
uses  a  seasonally  averaged  profile  (surface  to  bottom).  Although  the  applications 
if  salinity  differ,  tests  were  conducted  in  an  attempt  to  address  the  effects  on  the 
values  of  sound  speed  caused  by  the  two  different  approaches.  The  effects  on  sound 
speed  comparing  constant  salinity  with  variable  and  average  salinity  compared  with 
35.00  PPT  salinity  are  presented  in  tables  (  I  )  and  (II)  respectively. 


F.  Sound  Speed  Equation. A 

Common  to  both  the  ST  HAS  and  TCAPS  enviromental  data  output *  i*  a 
vertical  Aound  Apeed  profile.  The  value*  of  the  S THAS  *ound  Apeed*  profile* 
mau  differ  from  the  TCAPS  Aound  Apeed  profile*  becauAe  of  the  fcllov'iw ; 

Equation  uAed  to  compute  Aotmd  Aveed* 
ST  HAS  Lerotf'*  Fount  ion 

T CAP5  I''il*on'A  Equation 


( LeAOt /)  -  1.2  -  0. OUT-18)  ( equation  applicable  {ok  Aalinitie *  50-40  PET  and 

temperature*  0-30*). 

(WilAcn)  ^  =  1.398  +  3.384  x  10~3  (S-35) -1 . 12  x  10~2T  +  7.70  x  10 ~?P 

+  7.77  x  Tfl'V  -  1.29  x  T<fV  +  3.16  x  10~8PT 

-9  2 

+1.98  X  10  PT 

where  v  =  Aound  Apeed  (meter*/*ec) 

a  -  *  alinit.it  (  pt^t) 

T  =■  temperature  (<’C) 


Table  t  T  U )  Ia  taken  ( from  "An  Examination  of  the  Sen*itivitie*  of  A ccuAtic 
Performance  Prediction  (APP)  High  Level  Product*  to  Variation*  in  the  Ha  nine 
Environment" ,  bu  P.  Etter  and  V.  R.  F ruin,  MAP  Tccbucal  Report  *143,  April 
1980,  (SECRET).  In  the  opinion  of  the  author,  ba*ed  on  previou*  tooth  performed 
bu  Hr.  Paul  Etter  (MAP  Incorporated) ,  the  difference*  in  the  value*  oi  *cund 
Apeed  due  to  the  difference*  betioeen  Lerru'A  and  b'ilAon'A  equation*  are  veto 
Amall  and  not  Aignif leant. 
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G.  Sofitwivie.  Coding  E-Viosu  oh.  P-lobiem  Asizcu> 

During  their  supportive  investigation,  ODSI  had  noticed  some  programming 
errors  in  the  ICAPS  and  SIMAS  computer  codings.  The  specific  errors  uncovered  for 
ICAPS  and  SIMAS  have  been  documented  by  ODSI  and  are  described  in  Appendix  (C)  and  (D). 
Presented  below  from  each  system  are  a  few  selected  program  errors  which  are  considered 
noteworthy: 

Significant  SIMAS  Coding  Errors: 

1)  Metric  unit  inputted  data  are  converted  at  least  twice  into 
English  units  (i.e.,  5  meters/second  converted  to  16.4  feet/second 
then  reconverted  (16.4  m/s  x  3.3)  to  52.6  feet/second). 

2)  Erroneous  determination  of  the  layer  depth  by  execution  of  a 
DO  Loop  causing  erroneous  results,  see  appendix  (C). 

Significant  ICAPS  Coding  Errors: 

1)  When  the  use  of  metric  units  are  not  specified,  the  historical  bathymetry 
in  meters  in  incorrectly  multiplied  by  feet-to-meter  conversion  factor, 
see  appendix  (D). 

Program  "cautions"  where  potential  errors  can  occur  are  contained  on  Appendix 
(E)  and  (F). 

Coding  errors  described  in  Appendix  C-F  have  been  pointed  out  and  presented  to 
SIMAS  (NUSC)  and  ICAPS  (N00)  as  this  study  was  still  taking  place.  The  errors  may 
have  already  been  corrected  by  each  system  and  if  so,  did  not  affect  the  results  of 
this  study. 
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H.  Causeis)  o f  "4-£q>w.|^eau£"  dLmeAzn.au  in  the  SIMAS  and  ICAPS  sound  speed 
profiles  found  In  WPRflA  Technical  Note  *66. 

I n  an  attempt  to  determine  thz  cause  of  thz  "significant"  difference  found 
in  thz  sound  speed  profiles  of.  MORV A  TM  Mo.  66  titled,  "Comparison  of  the. 

ICAPS  and  SIMAS  Historical  Environmental  Data  Beau" ,  bu  F.  Hashimoto ,  dated 
28  September  7979,  there  were  eight  (8)  sites  randomlu  selected  from  that  studu 
for  sampling .  Each  of  thz  eiaht  sites  were  evaluated.  Evaluation s  bit  J. 

Look-in,  8.  Scaife  [Ocean  Data  Sustems,  Tnc. )  were  conducted  to  determine 
the  degree  of  influence  cauaed  bu  the  merge  algorithm,  the  historical  profile, 
the  bathuthermograph  data  on  bu  the  salinitu.  Table  [111]  bel.o w  represents 
a  matrix  showing  the  relative  signi  f  leant/ u  of  each  facto  A  at  each  site  location." 


Site 

Merge 

Algorithm 

Historical 

Profile 

Bathuthermograph 

Data 

Salinitu 

L'ariatlore 

Loccticn 

1C* 

Dominant 

Slight 

33. OX 

7  77. 5®r-' 

IF* 

Dominant. 

Slight 

47.  OX' 

743.  fow 

IF 

Dominant. 

Slight 

Slight 

<same  as 

above 

2  A* 

Dominant 

Moderate 

7  3.  OX 

38.0°(” 

2F 

Moderate 

Dominant 

Sliaht 

4  7.5X 

77, 5°f" 

2F 

Moderate 

Dominant 

Moderate 

same  as 

above 

2H* 

Dominant 

Moderate 

39.  FX 

68. 

5A* 

Dominant 

37.  CX 

17. F°F 

TABLE  ITT  OUALITATIVE  SIGNIFICANT  CAUSE  OF  DIFFIDENCES 
TM  SIMAS  I CAPS  COMPAPTSOMS  (OUST} 


*  In  these  SIMAS  generated  results,  the  BT  data  is  considered  as  probable  errered 
data,  so  that  the  generated  profile  is  the  historical,  profile  with  a  lauer  depth 
revision. 


V.  DISCUSSION 


If  one  were  to  glance  at  both  the  ICAPS  and  SIMAS  systems,  one  could 
conceivably  conclude  that  both  are  extremely  similar  and  do  not  differ  signi¬ 
ficantly.  This  conclusion  can  be  reached  when  the  following  are  lightly 
considered: 

•  Both  ICAPS  and  SIMAS  have  an  environmental  data  base 

•  Both  ICAPS  and  SIMAS  utilize  In-situ  BT  data  as  input 

•  Both  compute  a  final  sound  speed  profile 

•  Both  are  designed  for  "on-board"  and  "on-scene"  prediction 
systems 

Although  both  systems  have  and  do  all  of  the  above,  closer  investigations, 
as  presented  in  section  IV  and  in  the  Appendixes,  show  that  significant  basic 
differences  do  exist  and  strongly  indicate  that  the  ICAPS  and  SIMAS  systems 
are,  indeed,  not  identical  and  do  differ  markedly  in  some  areas. 

The  major  areas  where  each  system  differ  markedly  are: 

•  in  their  philosophies,  as  described  in  section  IV  (A) 

•  in  their  handling  of  the  In-situ  BT  data  INPUT  profile,  as 
described  in  section  IV  (B) 

e  in  the  historical  environmental  data  bases,  as  described  in 
section  IV  (C) 

•  -n  their  merge  methodologies  of  BT  to  the  historical  profile, 
as  described  in  section  IV  (D). 

The  areas  where  differences  exist  but  do  not  appear  to  be  significant  are: 

•  in  the  application  of  salinity  (between  approximately  34. 50  PPT 
to  36.00  PPT  above  1000  feet). 

•  in  the  sound  speed  equations 

The  areas  where  immediate  improvements  can  be  made  are. 

•  in  bathymetry,  through  the  use  of  a  larger,  more  recent,  hioher 
resolution,  automated  data  set  of  bathymetry.  See  Appendix  (G) 
for  description  and  details  regarding  resolution  and  coverage  of 
SYNBAPS  by  Mr.  R.  VanVfyckhouse  (N0RDA). 

•  in  the  historical  temperature,  salinity  and  sound  speed  profiles 
through  use  of  a  larger  (XBT,  AXBT  and  BT)  more-  recent  (up  to  1979), 
vertically  stable,  monthly  temperature,  seasonal  salinity  and 
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monthly  sound  speed  data  fields.  Excerpts  from  the  Pacific 
Sierra  Research,  Corp.  (PSR)  Report  £922,  "Evaluation  of 
Standard  Ocean  Candidates",  by  J.  Colborn,  S.  C.  Daubin,  Jr., 
E.  Hashimoto,  and  F.  J.  Ryan,  March  1980,  describinp  the 
Generalized  Digital  Environmental  Model  by  T.  Davis  (N00) 
are  presented  in  Appendix  (H). 
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V.  DISCUSSION 


If  one  were  to  glance  at  both  the  ICAPS  and  SIMAS  systems,  one  could 
conceivably  conclude  that  both  are  extremely  similar  and  do  not  differ  signi¬ 
ficantly.  This  conclusion  can  be  reached  when  the  following  are  lightly 
considered: 

•  Both  ICAPS  and  SIMAS  have  an  environmental  data  base 

•  Both  ICAPS  and  SIMAS  utilize  In-situ  BT  data  as  input 

•  Both  compute  a  final  sound  speed  profile 

•  Both  are  designed  for  "on-board"  and  "on-scene"  prediction 
systems 

Although  both  systems  have  and  do  all  of  the  above,  closer  investigations, 
as  presented  in  section  IV  and  in  the  Appendixes,  show  that  significant  basic 
differences  do  exist  and  strongly  indicate  that  the  ICAPS  and  SIMAS  systems 
are,  indeed,  not  identical  and  do  differ  markedly  in  some  areas. 

The  ma.ior  areas  where  each  system  differ  markedly  are: 

•  in  their  philosophies,  as  described  in  section  IV  (A) 

•  in  their  handling  of  the  In-situ  BT  data  INPUT  profile,  as 
described  in  section  IV  (B) 

•  in  the  historical  environmental  data  bases,  as  described  in 
section  IV  (C) 

•  in  their  merge  methodologies  of  BT  to  the  historical  profile, 
as  described  in  section  IV  (D). 

The  areas  where  differences  exist  but  do  not  appear  to  be  significant  are: 

•  in  the  application  of  salinity  (between  approximately  34.50  ppt 
to  36.00  PPT  above  1000  feet). 

t  in  the  sound  speed  equations 

The  areas  where  immediate  improvements  can  be  made  are: 

•  in  bathymetry,  through  the  use  of  a  larger,  more  recent,  higher 
resolution,  automated  data  set  of  bathymetry.  See  Appendix  (G) 
for  description  and  details  regarding  resolution  and  coveraqe  of 
SYNBAPS  by  Mr.  R.  VanWyckhouse  (NORDA). 

•  in  the  historical  temperature,  salinity  and  sound  speed  profiles 
through  use  of  a  larger  (XBT,  AXBT  and  BT)  more  recent  (up  to  1979), 
vertically  stable,  monthly  temperature,  seasonal  salinity  and 


monthly  sound  speed  data  fields.  Excerpts  from  the  Pacific 
Sierra  Research,  Corp.  (PSR)  Report  #922,  "Evaluation  of 
Standard  Ocean  Candidates",  by  J.  Colborn,  S.  C.  Daubin,  Jr., 
E.  Hashimoto,  and  F.  J.  Ryan,  March  1980,  describing  the 
Generalized  Digital  Environmental  Model  by  T.  Davis  (NOO) 
are  presented  in  Appendix  (H). 
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AND  DEVELOPMENT  ACTTVITY  NSTL  STATION  MS  E  HASHIMOTO 
30  SEP  80  N0R0A-TN-69 


UNCLASSIFIED 


F/G  9/2 


microcopy  resolution  test 

national  burlau  or  standards  1 


SIMAS(S) 

ICAPS(I) 

Constant 

Salinity 

Variable 

Salinity 

Sound 

Sound 

Site 

Depth 

(ft) 

Salinitv 

(PPT) 

Speed 

(ft/sec) 

Salinity 

(PPT) 

Speed 

(ft/sec) 

as 

(1-S) 

2.V 

(1-S) 

1C 

34.43 

4964.8 

34.41 

4964.9 

-0.02 

0.1 

S3  0. 1 

494S.5 

34.33 

494S.3 

-0.10 

-0.2 

1361. S 

4911.4 

34.07 

4910.1 

-0.36 

-1.3 

IF 
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33.40 
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-0.7 

767  . S 
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4887.6 

0.32 
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1630.7 

4852.9 

33.97 

4S54.4 

.  . 

0.31 

1.5 

2  A 

0.0 

36.34 

36.54 
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1.0 

734.9 

36.19 

4948.2 

-0.4 

1492.9 

4912.8 

35.33 

4908.9 

-1.01 

-3.9 

2H 

0.0 

33.72 

4917.2 

32.4  1 

4912.2 

-1.31 

-5.0 

807.1 

4920.3 

34.76 

4924.8 

1.04 

4.5 

1656.9 

4873.4 

34.91 

4878. S 

1.19 

5.4 

TABLE  I:  CONSTANT  VERSUS  VARIABLE  SALINITY 
EFFECTS  ON  SOUND  SPEED  PROFILES  (PDSI) 
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O  d 


UlL  <  I 


Site 


Salinity 

**500' 

(PPT) 

SIMAS  Data 

Base  Salinity(S) 
(PPT) 

AS 

S5Q0~S 

(PPT) 

Max  AV  est. 

(ft/sec) 

34.43 

35.0 

-.57 

-2.62 

33.66 

35.0 

-1.34 

-6.16 

36.34 

35.0 

1.34 

6.16 

33.72 

35.0 

-1.28 

-5.89 

Fcnn  =  Value  of  ICAPS  Salinity  at 
500  500  feet. 


TABLE  HjAVERAGE  VERSUS  35  PPT 
SALINITY  EFFECTS  ON  SOUND  SPEED  (ODSI) 
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Table  IV.  The  Maximum  Difference  ia  Sound  Speed  Among  Three  Frequently  Used  Sound 
Speed  Equations  (U)  ( T  MKi  *X  t  .»  » »  M/l  Cl  Time 


Point 

Depth 

Sound  Velocity 

(F/S) 

A  SVP 

* 

(ft) 

Leroy 

Frye  4  Pugh 

Wilson  (Rev) 

(F/S) 

1 

0 

5007.434 

5006 . 395 

5007.227 

1.039 

t 
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5009. 07S 

5008.082 

5008.871 

.996 

3 

200 
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5009.770 

5010.520 

.949 

4 

300 
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5010.969 
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.914 

5 

400 
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5012.168 

5012.844 

.875 
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5012 . "38 

5011. 8S7 

5012.545 

.8  51 

bOO 

5012.422 

5011.590 

5012.230 

.  S32 

S 

700 
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5011.906 

.813 

9 
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5010.957 

5011. 566 

.  797 

10 

900 

5011.598 

5010.617 
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.-81 

11 

1000 
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5010.266 
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.  761 

12 

1100 

5009. "46 
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.875 

13 

1200 

5008.953 
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5008.781 

.957 

14 

1300 

5008.664 

5007.652 

5008.492 

1.012 

15 

1400 

5007.848 

5006.801 

5007.676 

1.047 

16 

1500 

5007.539 

5006.477 

5007.367 

1.0b2 

1" 

1600 

5005.246 

5004. 1S4 

5005.078 

1 . 062 

18 

1800 

5002.414 

S 001 . 430 

5002.258 

.984 

19 

2000 

4998. "9" 

499". 867 

4998. b45 

.930 

:o 

2500 

4987.  t>29 

4986 . 863 

4987.508 

.  “66 

:i 

5000 

4966.602 

4966.023 

4966.543 

.  2"9 

■>  -i 

3500 

4933.992 

4953."34 

4933.996 

.  262 

25 

4000 

4911. "27 

4911. "15 

4911.-50 

.  015 

24 

4250 

4907.130 

4907.211 

4907.184 

.031 

25 

4500 

4903.258 

4905.320 

4903.254 

.  066 

26 

4750 

4901.484 

4901.562 

4901.480 

.082 

-)  “T 

5000 

4901.926 

4902.020 

4901.922 

.  098 

28 

5250 

4902.355 

4902.457 

4902.548 

.  109 

29 

5500 

4902.770 

4902.875 

4902.762 

.115 

50 

5750 

4905.465 

4905.578 

4905.461 

.  11" 

31 

6000 

4908. 164 

4908.281 

4908.160 

.  121 

32 

6500 

4915.109 

4915.234 

4915.109 

.  125 

33 

7000 

4922.074 

4922.215 

4922. 073 

.  141 

34 

8000 

4936.840 

4937.016 

4936.855 

.  1"6 

35 

9000 

4952.473 

4952.695 

4952.496 

■>22 

36 

10000 

4966 . 656 

4966.957 

4966. b95 

.  2S1 

37 

11000 

4980.941 

4981.289 

4980.992 

.  348 

38 

12000 

4996.094 

4996.516 

4996. 160 

.422 

39 

15000 

5011.352 

5011.848 

5011.430 

.496 

40 

14500 

5034.816 

5035.434 

5034.902 

.618 

94 

Mean 

.568 

APPENDIX  (A) 
SIMAS  LOGIC  FLOW 
by 

J.  Lock! in 
and 

B.  Scaiffe 
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Is  BT  data 

to  be  used? 

Yes 

No 

Input  BT  data  (depth, 
temperature)  pairs: 
(ZBi ,  TB^ ) ,  1=1, NB  «  24 


Set  SSP  being  genera¬ 
ted  to  be  the  histori¬ 
cal  profile. 


Are  any  points  to  be 
edited? 


Input  data  to  add,  de¬ 
lete  or  change  any 
point. 


METRIC  input  units. 


Convert  BT  points  to 
English  units. 


Via  LEROY's  equation, 
calculate  sound  speeds 
(VH.j)  using  BT  data 
and  historical  salini¬ 
ty.  VH.=  ^  (ZBi,TBi,S), 


Calculate  the  layer 
depth. 


Insert  the  souno  speed 
point  at  1000  ft.  into 
the  BT  SSP. 


B 
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m 
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APPENDIX  (B) 
ICAPS  LOGIC  FLOW 
by 

J.  Locklin 
and 

B.  Scaiffe 


I CAPS  Sound  Speed  Profile  (SSP)  Generation  (UNI VAC  Version) 


»)  »rci 
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Is  any  BT  profile 
depth  a  the 
bathymetry? 


Is  the  maximum  BT 
depth  greater  than  50 
meters  above  the 
maximum  historical 
profi le  depth? 


Find  index  I  to  ter¬ 
minate  BT  profile  at 
the  shallowest  BT 
depth  3  bathymetry. 
Set  N3T=I . 


Define  estimated  pro¬ 
file  (depth,  tempera¬ 
ture)  to  be  truncated 
BT  profile. 


Define  the  estimated 
profile  (depth,  tem¬ 
perature)  to  be  the 
BT  profi le. 


Select  the  index,  k, 
for  the  shallowest 
depth,  Zk,  in  the 
historical  profile 
where  the  maximum 
BT  depth  is  « 

Zk -50  meters. 


Initialize  the  first 
NBT  number  of  genera¬ 
ted  (depth,  tempera¬ 
ture)  profile  points 
to  be  the  BT  profi le. 


The  estimated  profile 
has  N=NBT  pairs  of 
points  (ZM.j,  TM.j) 
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Update  the  estimated  profile  so  that  the  deepest  depth  is 
the  bathymetry. 


Is  bathymetry  greater 
than  ZM^? 


Linearly  extrapolate 
temperature  and 
salinity,  and  add  a 
point  to  the  profile 
at  ZM^'bathymetry. 

Set  N*N+1. 


Linearly  interpolate 
temperature  and 
salinity,  and  set 
Zf^|=bathymetry. 


Correct  the  estimated  salinity  values  for  stability. 


m 


- * _ 

Loop  over  i=l,...,N. 
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APPENDIX  (C) 


2.  SIMAS 

2.1  PROGRAM  ERRORS 

2.1.1  Conversion  errors  occur  involving  subroutines  BT  and  EDITBT.  Metric 
inputs  are  converted  to  English  units  at  least  twice.  We  understand 
that  this  is  corrected  at  NUSC  but  not  at  NORDA. 

2.1.2  Subroutine  CKBT  shows  an  error  in  the  algorithm  for  modifying  the 
historical  data  to  reflect  the  BT  profile  layer  depth.  The  DO  loop 
at  FORTRAN  statement  number  550  may  be  executed  with  K=N15+],  i.e. 

DO  600  J=N15+1  ,N15 

thus  causing  erroneoud  results. 
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APPENDIX  (D) 

I.  ICAPS 

1.1  PROGRAM  ERRORS 

1.1.1  Univac  Version 

1.1. 1.1  The  synthesized  profile  contains  salinities  that  are  initialized  by 

interpolating  values  from  the  historical  profile,  and  when  a  depth  of 
the  synthesized  profile  is  equal  to  or  exceeds  the  maximum  depth  of 
the  historical  profile  (a  most  frequent  occurrence)  the  salinity  is  set  to 
be  the  historical  salinity  for  the  previous  depth.  This  can  be  corrected 
by  changing  the  inequality  ".LT."  to  ".LE."  in  the  IF  statement 
following  FORTRAN  statement  number  50  in  subroutine  MERGE.  This 
code  error  will  have  negligable  impact  for  the  deeper  profiles.  For 
shallower  profiles  the  impact  will  depend  upon  the  effect  on  velocity 
by  the  error  in  salinity. 

1.1. 1.2  When  the  user  wishes  to  use  only  the  historical  data  (no  BT  data: 
NDP=0  and  BOTZ  1.0)  and  metric  is  not  specified  (MOE  k  "M"),  the 
program  incorrectly  multiplies  the  historical  bathymetry  in  meters  by 
the  feet-to-meters  conversion  factor,  FTMT.  This  can  be  corrected  by 
inserting  after  format  statement  number  9001  in  the  main  program 

IF  (MOE  .NE.  1  HM)  BOTZ  =  FTMT*BOTZ 

and  deleting  the  statement  immediately  after  FORTRAN  statement 
number  10,  and  replacing  FORTRAN  statement  number  100  with 

100  CONTINUE. 

1.1. 1.3  An  error  message  is  not  utilized  when  appropriate  and  an  incorrect 
message  is  printed.  Change  the  line  just  preceeding  format  statement 
number  1045  to  read 

WRITE  (6,  1045)  IBTYP2. 

1.1. 1.4  In  the  main  program,  there  is  an  "MP"  variable  occurring  just  after 
format  statement  number  1013  that  is  not  initialized.  Apparently,  this 
was  to  have  been  defined  by  "NP"  in  the  preceeding  call  to  subroutine 
COMPUX.  This  is  corrected  by  changing  "NP"  to  "MP,"  the  last 
argument  in  the  argument  list  for  the  "CALL  COMPUX  ..."  statement, 
two  lines  after  FORTRAN  statement  number  360. 

1. 1.1.5  Subroutine  MERGE  calls  subroutine  SGMTST  to  perform  the  water 
column  stability  algorithm  from  the  bottom  of  the  merged  profile. 
Though  this  is  in  contrast  with  the  NOVA  version  that  starts  from  the 
bottom  of  the  BT,  it  is  more  reasonable  since  the  algorithm  is  applied 
just  below  the  BT  where  merging  occurs  (refer  to  paragraph  1.1 .2.2). 
Subroutine  SGMTST  has  an  incorrect  DIMENSION  statement  allocating 
insufficient  memory,  particularly  for  the  density  array.  The  current 
dimension  for  30  elements  is  not  consistent  with  the  calling  routine. 
The  arrays  should  be  dimensioned  for  50  elements  each,  since  the 
calling  routine,  MERGE,  identifies  the  arrays  as  for  the  merged 
profile.  This  is  corrected  by  changing  the  dimension  statement  in 
subroutine  SGMTST  to  read 

DIMENSION  T( 50),  S(50),  SIGMA(50) 
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1.1. 1.6 

1.1. 1.6.1 

1.1. 1.6.2 

1.1.2 
1. 1.2.1 
1. 1.2.2 


Under  normal  operating  system  usage,  the  density  and  stability  calcu¬ 
lations  will  be  incorrect.  Mixed  expressions  contain  integers  raised  to 
negative  integer  powers  which  result  in  zero  values. 

Subroutine  S1GMAT  has  an  equation  "B  =  ..."  that  contains  the  integer 
10  raised  to  a  negative  integer  power. 

Subroutine  SGMTST  has  the  equation  "DS  =  ..."  that  contains  two 
integer  10's  raised  to  negative  powers. 

Nova  S00  version 

Same  error  as  described  in  paragraph  1.1. 1.1  above. 

The  program  should  perform  the  water  column  stability  algorithm 
from  the  bottom  of  the  merged  profile  as  the  Univac  version  does 
(refer  to  paragraph  1.1. 1.5)  —  rather  than  from  the  bottom  fo  the  BT. 
Thus,  the  effects  on  density  (refer  to  paragraph  1.3. 4.2)  due  to 
temperature  revisions  by  merging  just  below  the  deepest  BT  depth 
would  be  accounted  for.  This  can  be  corrected  by  changing  in 
subroutine  MERGE,  the  FORTRAN  statement  just  after  statement 
number  100  to: 

N=NOPTM-l 


and  two  lines  later,  change  to  read 

K=NOPTM-3. 
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APPENDIX  (E) 


2.2  PROGRAM  CAUTIONS  (SIMAS) 

2.2.1  The  generated  SSP  arrays  are  allocated  for  50  points.  Subroutine 
INSERT  adds  points  to  the  generated  profile,  and  subroutine  ENVIRN 
extends  the  up  to  25  revised  BT  profile  points  with  historical  points. 
The  program  does  not  properly  check  on  this  limit.  Structurally,  depth 
array  element  2(51)  is  the  location  of  velocity  C(l),  the  surface 
velocity  which  could  be  changed  to  be  the  Z(51)  value. 

2.2.2  Subroutine  BT  can  input  up  to  25  BT  temperature  profile  points. 
Subroutine  CKBT  calls  subroutine  INSERT  to  augment  the  BT  points 
with  the  interpolated  temperature  point  at  1090  feet,  whenever  the 
input  BT  profile  does  not  have  a  1000  ft.  point.  Since  the  program 
allocates  space  for  25  BT  points,  an  overflow  can  occur,  so  that  the 
surface  temperature,  T(l),  value  is  replaced  by  the  overflow  depth 
value,  D(26). 

2.2.3  Whenever  the  BT  profile  is  augmented  by  a  point  at  1000  ft., 
extrapolation  is  performed  when  the  profile  is  to  be  extended  by  a 
1000  ft.  point.  The  SIMAS  program  also  does  not  consider  the 
reasonability  of  the  extrapolation,  see  paragraph  1.2. 1.2  above. 

2. 2. A  In  subroutines  ENVIRN  and  FORCST,  there  is  no  test  that  the 

corrected  bottom  depth  value  is  in  the  deepest  (last)  interval  which 
ends  with  the  input  fathometer  (or  chart)  value.  Consequently, 
extrapolation  may  occur. 

2.3  COMMENTS  —  The  methodology  is  to  shift  the  BT  SSP  so  that  at  1000 
fee:  the  shifted  BT  SSP  has  the  same  velocity  value  as  does  the 
historical  profile.  The  shifted  BT  SSP  is  then  extended  by  those 
historical  profile  points  below  the  deepest  BT  SSP  depth. 

.1  SIMAS  assumes  that  the  historical  profile  will  always  extend  below  the 
BT  profile.  This  assumption  can  only  be  valid  if,  over  each  geographi¬ 
cal  area,  the  associated  historical  profile  extends  down  to  the  deepest 
bathymetry  in  the  area,  and  any  deepest  BT  depth  value  would  not  be 
deeper  due  to  measurement  errors.  Though  the  assumption  will  usually 
be  valid,  it  is  too  strict  for  a  general  purpose  production  program. 
Since  the  program  does  not  verify  the  input  data  to  meet  this 
assumption,  errored  profiles  may  result  (refer  to  subroutine  ENVIRN 
lines  0052  through  0063), 

2.3.2  The  criteria  at  1000  feet  does  not  assure  a  smooth  transition  below 
that  depth  from  the  bottom  of  the  shifted  BT  profile  and  the  attached 
historical  points. 

2.3.3  Assuming  that  the  usage  of  the  layer  depth  is  important,  a  comment 
on  the  SIMAS  layer  depth  algorithm  is  appropriate.  Subroutine  LAVER 
defines  the  layer  depth  to  be  the  shallowest  depth  after  which  the 
velocity  gradient  first  becomes  negative.  This  algorithm  requires  that 
the  BT  data  not  be  noisy,  and  it  consists  of  points  that  only  define  the 
essential  shape  of  the  profile. 
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appendix  (F) 


1.2 

1.2.1 

1. 2.1.1 

1.2. 1.2 

1.3 

1.3.1 

1.3.2 

1.3.3 

1.3.4 

1.3.4. 1 


PROGRAM  CAUTIONS 
ICAPS,  Univac  and  Nova  versions 

When  a  merged  profile  is  created,  the  number  of  points  defining  the 
profile  is  not  verified  to  avoid  exceeding  the  allocated  storage  for  a 
maximum  of  50  points.  In  the  event  that  a  merged  profile  should  have 
more  than  50  points,  the  overflow  would  incorrectly  redefine  neces¬ 
sary  program  locations  (e.g.  the  near-surface  temperatures,  TM(l), 
TM(2),  ...,  would  be  redefined  by  depth  values  ZM(51),  ZM(52),  ..., 
respectively,  etc.) 

When  the  bathymetry  is  below  the  ICAPS  merged  profile,  the  merged 
profile  is  augmented  by  one  more  point  defined  by  the  bathymetry,  and 
temperature  and  salinity  values  extrapolated  from  the  deepest  (last) 
depth  interval  in  the  profile.  The  ICAPS  program  does  not  address  the 
question  "How  reasonable  are  the  extrapolated  values?".  Reasonabili¬ 
ty  of  these  values  is  dependent  upon  how  much  deeper  the  extrapola¬ 
tion  extends  the  profile  with  respect  to  the  depth  interval  being 
extrapolated  from,  and  the  gradient  on  that  same  interval. 

COMMENTS  —  ICAPS  is  a  merge  methodology  that  assumes  the 
bathythermograph,  BT,  data  is  the  best  current  information  upon 
which,  to  define  a  sound  speed  profile  (SSP);  and,  the  historical  data 
extending  below  the  BT  profile  is  the  most  reasonable  information 
available  to  merge  with  the  BT  data  to  form  representative  up-to-date 
SSP. 

A  shallow  water  case  is  defined  when  the  deepest  depth  of  the  input 
BT  data  is  above  200m.  This  excludes  a  merge  with  historical  data  — 
e.g.  an  operator  cannot  input  a  shallow  BT  profile  to  merge  with 
historical  data  for  a  deeper  water  mass.  In  shallow  areas  where  BT 
data  is  sufficient,  a  shallow  water  case  should  be  specified  by  another 
input  parameter,  rather  than  the  deepest  BT  profile  depth. 

For  a  shallow  water  case,  the  salinity  is  set  to  35  PPT.  This  is 
considered  valid  everywhere  except  in  the  Mediterranean  Sea  where  38 
PPT  is  more  appropriate.  In  keeping  with  the  philosophy  of  making  use 
of  the  "best"  information  available,  historical  data  (salinity)  should 
still  be  utilized  whenever  it's  available. 

Coding  of  constants  is  more  efficient  when  using  the  form 
+X.XXE+YY,  rather  than  +X.XX*10.**  (-YY)  (see  subroutines  SGMTST 
and  S1GMAT). 

There  are  inconsistencies  between  the  ICAPS  program  and  the  refer¬ 
ence  report,  "The  ICAPS  Water  Mass  History  File,"  by  Alvan  Fisher, 
Hr.,  May  1978. 

The  program  implementation  of  the  synthesized  temperature  al¬ 
gorithm  (see  Enclosure  5),  is  reasonable,  though  the  approximation  to 
the  synthesizing  algorithm  is  unnecessary. 
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1.3. 4. 2  The  density  (sigma-t)  stability  algorithm  described  in  Appendix  A 

(reference  (1))  does  not  correspond  to  the  implementation.  The 
document  states,  "...  salinity  inversion  must  coincide  with  a  tempera¬ 
ture  inversion  ...".  The  code  assumes  that  if  the  density  at  depth  Z.  is 
greater  than  the  density  at  the  depth  below,  Zj  j,  (starting  from  ihe 
bottom  of  the  profile),  then  recompute  the  salinity  at  2-  via  the 
stability  expression  and  use  this  "corrected"  salinity  value ^ with  the 
temperature  to  recalculate  the  density  at  depth  Z-.  There  is  no 
application  of  the  criteria  to  adjust  salinities  that^  are  "...  within 
temperature  inversions  that  are  more  than  the  temperature  maximum 
minus  0.25°C  at  the  lower  boundary  of  the  inversion  ..."  One  should 
note  that  the  stability  expression  was  derived  for  a  constant  tempera¬ 
ture  of  10°C  and  salinity  ranges  of  30  to  40  PPT.  This  algorithm  to 
"correct"  salinity  is  applied  whenever  the  density  does  not  monotonic- 
ally  increase  uith  depth,  and  over  the  temperature  ranges  of  the  BT. 
The  BT  temperatures  can  range  from  approximately  a°C  on  up  to 
27°C.  This  wide  range  of  temperatures  suggests  a  temperature 
varying  stability  algorithm  —e.g.  at  depth  Z-,  adjust  the  historical 
salinity  by  the  increment  indicated  -by  the  thange  in  tc-rnperature 
(merge  temperature  minus  historical  temperature)  to  stabilize  the 
water  column  at  depth  Z^ 

1.3.5  Both  program  versions  retrieve  a  historical  profile  as  described  in  the 
report  "Description  of  1CAPS  Environmental  Data  Structure"  by  ]ohn 
Lever,  NAVOCEANO  T.\  3700-S2-79.  March  1979.  The  report  presents 
the  file  structure  and  the  retrieval  algorithm.  The  following  com¬ 
ments  may  enable  the  report  to  be  more  useful:  particularly  for  profile 
revisions. 


1.3. 5.1  Expand  the  text  to  specifically  define  the  o-dering  of  profiles  (water 
masses)  for  a  geographical  location.  Although  the  ordering  is  intrinsic 
to  the  retrieval  design  (refer  to  report  figure  6),  the  following 
conditions  being  stated  would  clarify  the  profile  data: 

When  more  than  one  water  mass  represent  a  geographical  location,  let 
1<N€5  be  the  number  of  profiles,  then  XMlNT2(i),  X.MAXT2G),  i=i,  S' 
must  be  defined  (temperature  tolerances  at  200m),  and 

1)  XMINT2(i)<X.MAXT2(;)  for  i  =  l,  N 

2)  XMINT2G)  =  XMAXT2U-1)  for  i=2,  N 

must  be  true  for  continuity.  The  exception  to  this  is  that  two  adjacent 
profiles  can  have  the  same  tolerances,  X.MINT2(j)  =  XMINT2(j  +  l)  and 
XMAXT2(j)  =  XM  AXT2(j  + 1 ).  For  this  exception,  there  must  be  a  non¬ 
zero  flag  (NOGL(j))  in  order  to  consider  the  j+1  profile.  The  tie 
breaker  is  the  200-300m  gradient  tolerances  (XMINGL(j),  X.M AXGL(j)) 
and  (X MlNGL(j+ 1 ),  XMAXGL(j+l)).  When  NOGL(j)  is  non-zero,  then 
XMlNGL(j)  and  XMAXGL(j)  must  be  defined.  Profile  j  is  selected  for 
the  BT  temperature  gradient  less  than  XMAXGL(j),  otherwise  profile 
j+1  is  selected.  If  the  gradient  tolerances  are  given  for  j+1,  then 
XMAXGL(j)  =  XMINGL(j+  1)  must  occur  for  continuity.  In  practice,  the 
tolerances  at  j+1  are  only  necessary  to  make  an  output  message 
meaningful  when  the  BT  gradient  at  200-300m  is  greater  than 
XMAXGL(j  +  l).  B 

1.3. 5. 2  On  pages  3  and  9,  state  that  the  maximum  number  of  points  for  a 
historical  profile  is  45. 

1.3. 5. 3  On  pages  10  and  11,  the  references  labelled  "400"  should  be  labelled 
"80"  in  order  to  agree  with  the  programs,  as  do  the  other  labelled 
references. 

1.3. 5. 4  On  page  3,  line  2  of  the  second  paragraph,  "NOPTS(N)"  should  be 
"NOPTl(N)"  since  NOPTS  is  set  to  the  NOPTI(I)  value  for  the  selected 
profile. 

1.3.5  On  page  3,  the  text  should  state  that  NPRF^5. 
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APPENDIX  (G) 

SYNTHETIC  BATHYMETRIC  PROFILING  SYSTEM 
(SYNBAPS) 


BY 

Mr.  Roger  VanWyckhouse 

(Naval  Ocean  Research  and  Development  Activity) 
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Synthetic  Bathymetric  Profiling  System 

by 

Roger  J.  VanWyckhouse 
Ocean  Acoustics  Division 
Naval  Oceanographic  Laboratory 
Naval  Ocean  Research  and  Development  Activity 


The  Synthetic  Bathymetric  Profiling  System  (SYNBAPS)  is  a  combination 
cf  digital  computer  software  (programs)  and  a  random-access  storage  file 
:*  cricded  bathymetric  data,  employed  to  rapidly  generate  random,  great- 
c;-c:e,  bathymetric  profiles  suitable  for  acoustic  propagation  modeling. 

T  N5APS  is  comoletely  automatic,  requiring  only  the  input,  via  a  control 
car;,  cf  the  latitude  and  longitude  of  the  beginning  point,  the  bearing, 
ir.z  the  maximum  range  or  the  locations  of  the  beginning  and  end  points. 

Tre  generated  profile  is  available  in  two  forms.  The  first  is  a  computer 
t'-.-n  crofT.e  v>ere  range  in  whole  nautical  miles  is  plotted  against  depth, 
■'  eicoer  rezers,  'a thorns,  or  feet;  the  second  is  a  card  image  file  of  the 
:  fre  data.  The  prefile  outputs  in  the  card  image  file  are  available  on 
r-cr.rtic  case,  munched  cards  or  listing.  Single  depth  point  retrieval  is 
ar:ar.ev-  feature  of  the  system. 


•&  reed  for  a  computerized  bathymetric  data  bank  and  techniques  for 
'  mar- relating  large  quantities  of  data  became  evident  as  demand  for 
•or::  profiles  increased  and  became  more  urgent.  It  became  increasingly 
.To  to  satisfy  these  demands  through,  manual  compilation  of  depth 

contouring,  and  profile  constructions.  A  massive  recompilation 
iysis  of  bathymetric  data,  systematic  revision  cf  all  bathymetric 
:"ts  in  the  North  Atlantic  and  North  Pacific  Oceans,  including  extension 
.■erane  to  the  equator,  was  recently  completed.  The 

.•sing  the  existing  official  data  bank,  of  bathymetric 

i  or 


<  r_>  , 


-  ^  . 


u  cove! 

:icality  gt 

mgs  for  machine  generation  of  profiles  became  apparent.  The  need 
racial  iced  bathymetric  data  bank  to  support  acoustic  -  oceanographic 
■■a  gave 


■'-Mng  gave  rise  to  development  of  a  synthetic  bathymetric  profiling  project 


me 


)c  the  new  bathymetric  contour  charts  as  the  data  base.  The  project 
Toped  procedures  for  digitizing  the  contour  charts,  and  computer  programs 
subroutines  for  data  storage  and  retrieval  and  for  profile  generation. 


The  SYNBAPS  data  base  was  designed  to  meet  the  specific  and  immediate 
need  for  bathymetric  profiles  for  acoustic  modeling.  However,  properly 
used,  the  data  base  offers  a  myriad  of  applications  beyond  its  preliminary 
design.  The  coverage  of  the  data  base  is  shown  in  figure  1. 

Often  in  naval  planning  as  well  as  irt  naval  operations,  speed  is  as 
important  as  accuracy  when  information  is  needed.  SYNBAPS  is  not  ideally 
suited  to  hydrographic  charting  because  some  high  resolution  detail  is  lost, 
but  it  provides  very  rapid  responses.  SYNBAPS  has  these  additional  features: 
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•  Only  gridded  depth  points  are  stored  in  the  data  bank  (3,489,066  points), 

•  The  locations  of  the  depth  points  are  logically  structured  on  a  Mercator 
projection  by  5-longitudinal  min.utes/5-meridianal  part  intersections, 

•  Random  access  to  the  data  is  by  large  blocks  (691  five-degree  squares), 

9  The  data  bank  is  updated  by  replacing  the  blocks  of  data, 

•  The  size  of  the  data  bank  is  fixed  once  it  has  been  created  for  an 
ocean  area, 

«  Rapid  contouring  of  five-degree  squares  of  data  for  inspection  or 
editing  can  be  done  directly  for  the  Mercator  projection,  the  most 
commonly  used  projection  for  bathymetric  charts, 

t  Classified  survey  data,  in  chart  form,  can  be  incorporated  in  the  data 
base  with  no  compromise  of  security,  the  present  data  base  is 
unclassified, 

9  Highly  compacted  forms  of  the  accessing  program  and  the  data  base 
can  be  used  on  shipboard  or  other  platforms, 

9  Subsets  of  the  data  base  can  be  used  in  other  model  data  bases 
requiring  bathymetry  (a  derivative  of  SYNBAPS  with  1/6°  resolution 
is  being  used  in  the  Automated  Signal  Excess  Prediction  System 
( AS  EPS)  Data  Base), 

t  The  data  base  is  available  on  one  9  track,  1600  bpi,  EBCDIC  code 
magnetic  tape, 

e  The  data  base,  a  profile  accessing  program  and  an  editing  program 
are  fully  documented  and  provided  with  the  system, 

•  In  addition  to  the  documentation  provided  with  the  data  base  the 
following  references  are  available; 
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APPENDIX  (H) 

GENERALIZED  DIGITAL  ENVIRONMENTAL  MODEL 

(GDEfl) 

BY 

Dr.  T.  Davis  (Naval  Oceanographic  Office) 
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Overview 


"This  report  describes  and  evaluates  eight  existing  or  proposed 

oceanographic  models  as  candidates  for  Standard  Ocean,  a  data  retrieval 

system  to  be  installed  in  the  Lonq  Range  Acoustic  Propagation  Project 

(LRAPP)  data  bank.  The  primary  purpose  of  Standard  Ocean  is  to  provide 

range-dependent  sound-speed  profiles  for  input  to  NORDA's  numerical  acoustic 

models.  Standard  Ocean  will  also  be  used  to  support  the  objective  analysis 

of  environmental  data  collected  during  exercises  at  sea.  The  candidate 

systems  and  their  parent  organization  are  as  follows: 

AUTO-OCEAN  (NORDA)  GFOL  (Geophysical  Fluid  Dynamics 

Laboratory,  NCAA,  Princeton 
FIB/EOTS/EXTRA  (FNWC)  University) 

GDEM  (NAVOCEANO)  HYDAT  (FNWC) 

ICAPS  (NAVOCEANO)  SIMAS  (NUSC/New  London) 

ODSI  (Ocean  Data  Systems,  Inc.) 

The  Standard  Ocean  Evaluation  Group  assessed  each  candidate  according 
to  criteria  indicated  in  the  following  description  of  desired  Standard  Ocean 
capabilities.  Standard  Ocean  is  to  provide  accurate,  realistic,  and  seasonal 
(preferably  monthly)  surface-to-bottom  profiles  of  sound  speed,  temperature, 
ar.d  salinity  in  each  oceanic  1°  x  1°  square.  The  sound-speed  profiles  should 
be  in  a  format  suitable  for  numerical  acoustic  models.  The  profiles  should 
accurately  reproduce  all  acoustically  significant  features.  The  degree  of 
oceanic  variability  in  each  square  should  be  indicated.  Standard  Ocean  should 
operate  rapidly  and  inexpensively;  it  should  be  easily  usable  by  the  nonspecialist. 
Finally,  the  candidate  chosen  should  be  competitive  in  acquisition  cost  and 
avail  ability." 
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Generalized  Digital  Environmental  Model  (GDEH)  by  Dr.  T.  Davis,  NOO 
"GDEM  is  an  objective-analysis  program  that  produces  fields  of  analyzed 
temperature,  salinity,  and  sound  speed  on  a  30'x30'  orid  over  various  depth 
ranges.*  Preliminary  analysis  has  been  completed  for  the  North  Pacific  and 
Mediterranean  and  is  heginnino  for  the  North  Atlantic.  The  Indian  Ocean 
seoment  is  to  be  completed  sometime  in  the  future.  Observational  data  are 
insufficient  to  run  GDEM  in  the  southern  hemisphere. 

The  GDEM  component  for  the  upper  layer  models  seasonal  temperature  between 
the  surface  and  the  meroe  depth  (400m) ,+  annual  salinity  between  the  surface 
and  either  400m  in  the  Mediterranean  or  800m  in  the  Pacific,  and  seasonal 
sound  speed  from  the  surface  to  the  merge  depth  in  both  oceans.  The  sound 
speed  is  derived  from  the  analyzed  temperature,  salinity  fields. 

Below  the  surface  models  of  temperature,  salinity,  and  sound  speed, 
there  is  a  middepth,  two-season  (winter-summer)  sound-speed  model  from 
200m  to  2450m.  It  is  being  augmented  in  the  North  Atlantic  by  a  mid-depth 
temperature  and  salinity  model. 

The  lowest  component  is  an  annual  sound-speed  model  extending  from 
2000m  to  the  bottom.  There  is  no  evidence  that  the  developer  intends  to 
replace  this  model  with  a  temperature  or  salinity  model  in  the  near  future. 

Of  all  Standard  Ocean  candidates  using  objective-analysis  techniques, 

GDEM  comes  closest  to  reproducing  the  significant  features  in  vertical  sound- 
speed  profiles.  Comparison  of  model  profiles  v'ith  typical  observed  profiles 


♦This  model  has  received  our  most  detailed  scrutiny  to  date. 

+Merge  depth,  a  variable  peculiar  to  each  ocean  basin,  is  the 
depth  at  which  two  adjacent  layers  of  the  model  are  merged.  In  the 
early  1979  version  of  the  model  for  the  North  Pacific,  for  example, 
that  depth  was  fixed  for  the  entire  basin  at  400m.  We  understand 
that  the  latest  versions  of  the  model  adjust  merge  depth  over  the 
basin  to  fit  the  observed  data. 
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Model 

GDEM 


for  the  North  Pacific  showed  that  GOEM  maintains  horizontal  continuity  yet 
preserves  frontal  structure  in  some  cases.  It  appears  to  work  well  in  areas 
with  sparse  data  except  the  southern  hemisphere.  The  30’x30'  qrid  snacinq 
is  the  smallest  of  any  model  yet  evaluated.  Operation  is  rapid  (3  min  to 
edit  an  entire  tape  for  one  season  in  the  North  Pacific)  and  reported  to  be 
simple.* 


Temperature 


Salinity 


Sound  Speed _ Coverage 


Seasonal,  30'x30',  surface 
to  400m  at  standard  depths 

monthly  also  available 


Annual,  30'x30',  sur¬ 
face  to  800m  at  stan¬ 
dard  depths 

seasonal  also 
available 


Spatial  resolution, 
30'x30'.  Temporal 
resolution  (3  models): 

(1)  seasonal,  0-400m, 

(2)  2  seasons,  200- 
2450m,  (3)  annual, 
2000  m- bottom 


Medi ter  - 
ranean,  N 
Pacific, 
North 
Atlantic 


monthly  also  available 
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